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Abstract 

The HI Collaboration at the DESY electron-proton collider HERA has observed, 
in photoproduction and neutral-current deep-inelastic scattering, an unexpected 
excess of events with isolated leptons and missing transverse momentum, especially 
at large values of hadronic transverse momentum — a signature typical for single 
top-quark production. This observation is being substantiated in the HERA II 
run. Motivated by this, we evaluate the cross section of single top-quark photo- 
and electroproduction both in the standard model and its minimal supersymmetric 
extension, considering both minimal and non-minimal flavour-violation scenarios in 
the latter case. 

PACS numbers: 12.60.Jv, 13.60.Hb, 13.85.Ni, 14.65.Ha 



*Present address: Deutsches Elektronen-Synchrotron (DESY), Notkcstr. 85, 22607 Hamburg, Ger- 
many. 



1 Introduction 



Searches for single top-quark production via the neutral current (NC), 

e ± p^e ± t + X, (1) 

have been performed by the HI Collaboration [1] and the ZEUS Collaboration [2] at 
the DESY electron-proton collider HERA. The HI Collaboration found several events, 
leading to a cross section of a = 0.29^°'^ pb. Alternatively, assuming that the observed 
events are due to a statistical fluctuation, upper limits of 0.55 pb on a and of 0.27 on the 
anomalous twy coupling K tui were established at the 95% confidence level (CL). On the 
other hand, the ZEUS Collaboration found no evidence for top-quark production and was 
able to place upper bounds of 0.225 pb on a and of 0.174 on K tul at 95% CL. 

Single top-quark production via the charged current (CC) is possible at the tree level, 
but its cross section is less than 1 fb [3]. Furthermore, such events can be separated 
experimentally due to the absence of a scattered electron or positron in the final state. 

In this work, we calculate the cross section of process ([1]) in the standard model (SM) 
and its minimal supersymmetric (SUSY) extension (MSSM), considering both scenarios 
with minimal flavour violation (MFV) and non- minimal flavour violation (NMFV). The 
bulk of the cross section is due to photoproduction and electromagnetic deep-inelastic 
scattering (DIS), while the contribution due to the exchange of a virtual Z boson is 
greatly suppressed by its mass. Since photonic interactions cannot change flavour and 
the top quark does not appear as a parton in the proton, we are dealing here with a 
loop-induced process. In the SM and the MFV MSSM, its cross section is further sup- 
pressed by the smallness of the contributing elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. In more general MSSM scenarios, misalignment between the quark and 
squark sectors can appear, and the CKM matrix is no longer the only source of flavour 
violation. Thus, the flavour-changing (FC) couplings are not Cabibbo suppressed, and 
sizeable contributions to FC NC processes can occur. On the other hand, there are strong 
experimental bounds on squark mixing involving the first generation, coming from data 
on K°—K and D°-Tf mixing [I]. If the squark mixing involving the first generation is 
neglected completely, there are no new vertices for incoming up-quarks compared to the 
MFV MSSM. Thus, every new contribution in this scenario is suppressed by the parton 
distribution function (PDF) of the charm quark in the proton. Because the SM cross 
section for NC single top-quark production is highly suppressed, every detected event is 
an indication of physics beyond the SM. 

This paper is organised as follows. In Section we describe the analytical calculation 
of the SM cross section. In Section [3l we outline the theoretical framework of FC inter- 
actions in the MSSM. The numerical analysis is presented in Section HI Our conclusions 
are summarised in Section [51 

2 SM cross section 

In this section, the analytical calculation of the SM cross section is described. As indicated 
in Fig.[U we denote the four-momenta of the incoming proton, parton (up or charm quark), 
and electron by P, p, and k, respectively, and those of the outgoing top quark and electron 
by p' and k', respectively. As for the centre-of-mass (CM) energy and the top-quark mass, 
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we have S = (P + k) 2 and mf = {p') 2 . We neglect the masses of the proton, the incoming 
quarks, and the electron, so that P 2 = p 2 = k 2 = k' 2 = 0. The four-momentum of the 
exchanged photon is given by q = k — k', and, as usual, we introduce the virtuality variable 
Q 2 = —q 2 > 0. The variable y — (q ■ P)/(k ■ P) measures the relative electron energy loss 
in the proton rest frame. 

2.1 Electroproduction 

Single top-quark production in NC DIS occurs via the partonic subprocess 

e ± q -> e ± t, (2) 

where q = u,c. The Feynman diagrams contributing in the SM to process (J2]) with q = u 
are depicted in Fig. The ones for q = c are similar. The amplitude of this processes was 
also calculated in Ref. [5]. There appears at least one off-diagonal element of the CKM 
matrix in each term. If the contributions of the inner-quark flavours to a single Feynman 
diagram are added up and their masses are neglected, the amplitude of this Feynman 
diagram vanishes due to the unitarity of the CKM matrix. Thus, we cannot neglect the 
inner-quark masses here. Although our choice m c = < appears unphysical at first 
sight, it is inconsequential in practice. In fact, we verified that the use of a realistic value 
of m c affects our numerical results only insignificantly. A detailed proof of hard-scattering 
factorisation with the inclusion of heavy-quark masses may be found in Ref. [6] . 

Our calculation proceeds along the lines of Ref. [TJ. The differential cross section of 
process (T2]) reads 

where £, defined as p = £P, is the fraction of the proton momentum passed on to the 
incoming quark and \A4\ 2 is the squared amplitude averaged (summed) over the spin and 
colour degrees of freedom of the initial-state (final-state) particles. The kinematically 
allowed ranges of Q 2 and y are 

Qcut < Q < UmaxS — m t , 
Q 2 + m t 

^-^ <V< Z/max, (4) 

where Q 2 nt defines the demarcation between photoproduction and electroproduction and 
y max is an experimental acceptance cut. In Ref. PQ, values for Q 2 nt and y max are not 
specified. In our numerical analysis, we employ the typical values Q 2 nt = 4 GeV 2 and 
Umax — 0.95. Our numerical results are insensitive to the precise choice of y max as long it 
is close to unity. On the other hand, the dependence on Q 2 nt approximately cancels out 
in the combination of photoproduction and electroproduction, as we explicitly verify. 

In the parton model of QCD, the hadronic cross section of process (DQ) is obtained by 
convoluting the partonic cross section of process (T5]) with the appropriate PDF F q (^,fj,p), 
where \ip is the factorisation scale, and summing over q — u,c, as 



, d Q 2d yJhMr ' \ d Q 2d y/ pai t 
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There are two candidate mass scales for /j,p, namely \[Q 2 ~ and m t , and the optimal choice 
is likely to lie somewhere in between, at fip = {^fQ 2 + m^j /2 say. At any rate, we have 
fiF ^> m c , so that charm is an active quark flavour in the initial state, contributing at full 
strength via its PDF. 

The expression for \Mfof process (J2]) may be decomposed into a hadronic tensor H^" 
and a leptonic tensor L^, as 

-^2 e 2 



Q 

where e is the positron charge. We have 



\M\' = — 4 L^H^ (6) 



where 



H,u = \Y, H l H »i (7) 

spins 



Hp = u{p') (F^Pl + F 2lfl P R + F 3Pfl P L 

+ F 4Ptl P R + F 5 p'P L + F 6 p'P R ) u(p), 



with helicity projectors Pl,r = (1 =F 7 5 )/2 and form factors Fx, ... , Fq, which follow from 
the explicit evaluation of the Feynman diagrams shown in Fig. [2] and their counterparts 
for an incoming charm quark. The leptonic tensor may be decomposed into transverse 
and longitudinal components, as 



" ^ v — 

,2 



^{[l + (l-Z/) 2 ]^-4(l-y)ef}, (9) 



where 



AO 2 2 
(Q 2 + m 2 t ) Q 2 + mf 



^-M^f^) G?tV + ^' (10) 



To obtain the transversal and longitudinal parts of the cross section, the hadron tensor is 
contracted with the transversal and longitudinal parts of the lepton tensor, respectively. 

We generate and evaluate the Feynman diagrams in Fig. [21 with the virtual-photon leg 
amputated, with the help of the program packages FeynArts [8119] and FormCalc [9p0] . 
We work in 't Hooft-Feynman gauge and use dimensional regularisation to extract the 
ultraviolet (UV) divergences. We perform the Passarino-Veltman reduction and calculate 
the squared amplitude using the program package FeynCalc [11]. For the numerical 
evaluation of the standard scalar one-loop integrals, we employ the program package 
LoopTools [T0ffT2] . We perform the numerical integration with the aid of the program 
package Cuba [13]. As a check, we also calculate the amplitude of process with the 
help of the program package FeynCalc and by hand. All three independent calculations 
are found to lead to the same result. Furthermore, we verify the cancellation of the UV 
divergences, current conservation, and the reality of the squared amplitude. 
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2.2 Photoproduction 



In the photoproduction limit, the virtual photon is considered as real, with Q 2 = 0, so 
that the longitudinal part of the cross section vanishes. In turn, its energy distribution is 
described in the Weizsacker- Williams approximation by the electron-to-photon splitting 
function 



SM = £ 



l + (l-y)* Ql t , -/ 1 1 



In ^p- + 2ym, 



11) 



V Qmin \ Qcut Qmin , 

where a = e 2 /(47r) is Sommerfeld's fine-structure constant, m e is the electron mass, and 
Qmin = l/ 2m e/(l — y) corresponds to the kinematic lower bound. Thus, the cross section 
of process §2$ in photoproduction is given by 



da\ 

Jy) 



) = /r(vWv). (12) 

part 



with 

a, = ^\M,\ 2 5(ZyS-m 2 t ), (13) 
where M. 1 is the amplitude of —* t. The kinematically allowed range of y is 

2 

Try 

-J- < y < 2/max- (14) 

The hadronic cross section is again obtained by convoluting the partonic cross section for a 
given incoming quark with the corresponding PDF and summing over the incoming-quark 
flavours. Here, the only candidate mass scale for \ip is of order mt, and it is plausible 
to choose fMp = m t /2 so that there is a smooth transition between photoproduction and 
electroproduction. Again, we have fip 3> m c , so that the use of a charm PDF is justified. 



3 Minimal and non-minimal flavour violation in the 
MSSM 

In the MSSM, there are two sources of FC phenomena [14] • The first one is due to 
flavour mixing in the quark sector, just as in the SM. It is produced by the different 
flavour rotations in the up- and down-quark sectors, and its strength is driven by the off- 
diagonal CKM matrix elements. This mixing produces FC electroweak (EW) interaction 
terms involving CCs, now also involving charged Higgs bosons, and SUSY EW interaction 
terms of the chargino-quark-squark type. Thus, the SM Feynman diagrams of Fig. [2] are 
supplemented by those shown in Fig. [31 in which either charged Higgs bosons and down- 
type quarks or charginos and down-type squarks circulate in the loops. In the MFV 
MSSM, this is the only source of FC phenomena beyond the SM. 

The second source of FC phenomena, which is present in the NMFV MSSM, is due to 
the possible misalignment between the rotations that diagonalise the quark and squark 
sectors. When the squark mass matrix is expressed in the basis where the squark fields 
are parallel to the quark fields (the super CKM basis), it is in general non-diagonal 
in flavour space. This quark-squark misalignment produces new FC terms in NC as 
well as in CC interactions. In the SUSY QCD sector, the FC interaction terms involve 
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NCs of the gluino-quark-squark type. In the case of process (j5J), this gives rise to the 
additional Feynman diagrams shown in the first row of Fig. HI In the SUSY EW sector, 
the FC interaction terms involve NCs of the neutralino-quark-squark and the chargino- 
quark-squark type. The first type appears exclusively due to quark-squark misalignment, 
as in the SUSY-QCD case, whereas the second type receives contributions from both 
sources, quark-squark misalignment and CKM mixing. The additional Feynman diagrams 
involving neutralino-quark-squark interactions are displayed in the second row of Fig. HI 
In order to simplify our analysis in the NMFV MSSM, we take the CKM matrix to 
be diagonal, so that SM contributions (cf. Fig. [2]) and the genuine MFV-MSSM contri- 
butions, i.e. the charged-Higgs-quark contributions and the part of the chargino-squark 
contributions due to CKM mixing (cf. Fig. [3]), are zero. We are then left with the gluino- 
squark and neutralino-squark contributions (cf. Fig. @J and the residual chargino-squark 
contributions. 

Furthermore, we assume that the non-CKM squark mixing is significant only for tran- 
sitions between the second- and third-generation squarks, and that there is only left-left 
(LL) mixing, given by an ansatz similar as in Ref. [T5], where it is proportional to the 
product of the SUSY masses involved. This assumption is theoretically well motivated by 
the flavour-off-diagonal squark squared-mass entries that are radiatively induced via the 
evolution from high energies down to the EW scale according to the renormalisation group 
equations (RGEs) [16]. These RGEs predict that the FC LL entries scale with the square 
of the soft-SUSY-breaking masses, in contrast with the left-right (LR) or right-left (RL) 
and the right-right (RR) entries, which scale with one or zero powers, respectively. Thus, 
the hierarchy LL 3> LR, RL 3> RR is usually assumed. The same estimates also indicate 
that the LL entry for the mixing between the second- and third-generation squarks is the 
dominant one due to the larger quark-mass factors involved. On the other hand, the LR 
and RL entries are experimentally more constrained, mainly by b — * 57 data [T7]. With 
the previous assumption, the squark squared-mass matrices in the (v,L,CL,tL,UR,CR,tR) 
and (di, §l, b\, d,R, Sr, b R ) bases can be written as follows: 



( M lu 








m u X u 





\ 





Ml, c 


A' LL M ijC M Lit 





m c X c 








\i L M LiC M L , t 
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(15) 



(16) 



where 

M 2 La = M% q + m\ + cos(2/3)Ml (if - Q q s 2 w ) , 

M\ q = Ml q + m 2 q + cos(2f3)Q q s 2 w M 2 z (q = u, c, t), 

M\ q = Ml q + m 2 q + cos(2(3)Q q s 2 w M 2 z (q = d, s, b), 



6 



X q = A q - ^(tan/?)" 2 ^, 



(17) 



and \ l LL and X b LL measure the squark flavour mixing strengths in the t—c and b-s sectors, 
respectively. As for the SM parameters, m q , T 3 , and Q q are the mass, weak isospin, 
and electric charge of quark q; Mz is the Z-boson mass; and s w = sm6 w is the sine of 
the weak mixing angle 9 W . As for the MSSM parameters, tan/3 = v 2 /vi is the ratio of 
the vacuum expectation values of the two Higgs doublets; \i is the Higgs-higgsino mass 
parameter; A q are the trilinear Higgs-sfermion couplings; and Mq , My , and are 
the scalar masses. Owing to SU(2)l invariance, we have Mq u = Mg d , Mq- c = Mq s , and 
Mq t = Mn b . Further MSSM input parameters include the EW gaugino masses Mi and 
M 2 , the gluino mass M 3 , and the mass M40 of the CP-odd neutral Higgs boson A . 

In order to reduce the NMFV-MSSM parameter space, we make the following sim- 
plifying assumptions. We assume that the flavour mixing strengths in the t—c and b—s 
sectors coincide and put A = \ l LL = \ b LL for a simpler notation. Obviously, the choice 
A = represents the case of zero squark flavour mixing. We assume that the various tri- 
linear Higgs-sfermion couplings coincide and write Aq = A u = A c = A t = Ad = A s = At,. 
We assume that the scalar masses coincide thus defining the common SUSY mass scale 
M = Mq q = Mjj , c t x = Mf) { dsb y As for the gaugino masses, we impose the GUT rela- 
tion Mi = (5/3) (s 2 u /c 2 u )M 2 , where c 2 = 1 — s 2 , while we treat the gluino mass parameter 
M 3 as independent. We are thus left with eight independent MSSM parameters, namely, 
tan/3, M40, M , M 2 , M 3 , Aq, //, and A. 



4 Numerical analysis 

In this section, we present our numerical results. We adopt the SM parameters from 
Ref. [TH] and the effective masses of the down-type quarks from Ref. [T9] : 

a = 1/137.035 99911, a, = 0.1176, 

(he) 2 = 0.389 379 323 GeV 2 mb, m e = 0.510 998 918 MeV, 

M w = 80.403 GeV, M z = 91.1876 GeV, m t = 172.7 GeV, 

m d = 0.041 GeV, m s = 0.15 GeV, m b = 4.5 GeV. (18) 

Here, m t and m b correspond to pole masses, while and m s were determined so that their 
insertion in the perturbative formula for the vacuum polarisation function of the photon 
reproduces the result extracted from the total cross section of hadron production in e + e~ 
annihilation, measured as a function of the CM energy, using a subtracted dispersion 
relation. This set of quark masses is especially appropriate for quantitative studies in 
electroweak physics and is frequently employed in the literature. It is used here for 
definiteness and convenience. As a matter of principle, the precise definition of quark 
mass is not yet fixed in an analysis of leading order (LO) in QCD like ours. In fact, the 
freedom of choice of quark-mass definition contributes to the theoretical uncertainty. We 
employ the standard complex parametrisation of the CKM matrix in terms of three angles 
#12, 023, 013 and a phase 5 13 , 



V 



( ci 2 c 13 S12C13 s 13 e~ i5ia 

-S12C23 - c 12 s 2 3S 13 e i5ri c 12 c 23 - s 12 s 23 s 13 e t5is s 23 c 13 j. (19) 
V s 12 s 23 - c 12 c 23 s 13 e tSl3 -c 12 s 23 - s 12 c 23 s 13 e lSl3 c 23 c 13 
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where Cy = cos% and = sin%, and adopt from Ref. [18] the values 

si 2 = 0.2272, s 23 = 0.0422, s 13 = 0.0040, 5 13 = 1.00. (20) 

As for the proton PDFs, we employ the LO set CTEQ6L1 [20] by the Coordinated 
Theoretical-Experimental Project on QCD (CTEQ) Collaboration. We choose the fac- 
torisation scale to be fj,p = {^/Q 2 + m^j /2, with the understanding that Q 2 = in the 
case of photoproduction. At HERA II, electrons or positrons of energy E e = 27.6 GeV 
collide with protons of energy E p = 920 GeV in the laboratory frame, yielding a CM 
energy of = 319 GeV. 

4.1 Standard Model 

Figures [5] and [6] refer to the SM. Figure [5] shows the total cross section of process ([T|) as well 
as its photoproduction and electroproduction components as functions of y/S. We observe 
that, for our choice of Ql ut , the contribution due to photoproduction is approximately 
twice as large as the one due to electroproduction. At the CM energy of HERA, the 
cross section is of order 10~ 10 fb and thus many orders of magnitude too small to be 
measurable. A possible future electron-proton supercollider that uses the HERA proton 
beam with energy E p = 920 GeV and the electron beam of the international linear e + e~ 
collider (ILC) with energy E e = 500 GeV would have a CM energy of 1357 GeV. At this 
energy, the cross section is of order 10~ 8 fb and likewise not measurable. 

Figure [6] shows the Q 2 distribution of the electroproduction cross section as well as its 
transversal and longitudinal parts. We observe that the transversal part makes the major 
contribution. 

4.2 MSSM with minimal flavour violation 

For defmiteness, we assume that SUSY is broken according to the mSUGRA scenario of 
a Grand Unified Theory (GUT). We assign the following default values to the mSUGRA 
input parameters at the GUT scale: 

tan = 56, m = 1.25 TeV, m 1/2 = 140 GeV, 

A = -260 GeV, sign(» = +1, (21) 

where m is the universal scalar mass, mi/ 2 is the universal gaugino mass, and A is the 
universal trilinear Higgs-sfermion coupling. These values approximately maximise the 
cross section of process ([T]) and are in accordance with the experimental bounds from 
b — ► 57 decay and on the masses of the various SUSY particles [18]. We calculate the 
MSSM mass spectrum with the help of the program package SuSpect [21]. 

Figures I7T4T01 all show the total cross section of process ([T]) in the MFV MSSM. Specif- 
ically, Fig. [7] displays the y/S dependence, also separately for the photoproduction and 
electroproduction contributions, while Figs. [5HTU1 exhibit, for HERA experimental con- 
ditions, the dependencies on tan/?, mi/ 2 , and A , respectively, also separately for the 
charged- Higgs-boson and chargino contributions. 

From Fig. [7] we observe that the cross section is of order 10~ 5 fb at HERA energy 
and of order 10~ 3 fb for the future electron-proton supercollider mentioned above. Both 
values are too small to yield measurable results. 



8 



From Fig. [8] we learn that, as tan/? approaches its upper limit, the cross section 
strongly increases and is mainly generated by the loop diagrams involving charged Higgs 
bosons. For small values of tan /3, the SM contribution (w 10~ 10 fb) is dominant. Negative 
interference effects between the charged-Higgs, chargino, and SM contributions can be 
seen at large values of tan (3. 

From Fig. Owe see that, as approaches its lower limit, the cross section strongly 
increases and is essentially made up by the charged-Higgs-boson contribution alone. The 
latter is dominant throughout the whole mu^ range, but there are negative interference 
effects for all values of mu2- The significant suppression of the chargino contribution for 
tan j3 = 56 familiar from Fig. [8] is actually present for all values of m^. 

As is evident from Fig. [10j the cross section is largest for Aq = —260 GeV and falls off 
by one (two) orders of magnitude as Aq reaches —1 TeV (1 TeV). However, the variation 
with Aq is less significant than those with tan/? and mm. For tan/? = 56, the chargino 
contribution is several orders of magnitude smaller than the charged-Higgs one and almost 
independent of Aq. 

We conclude that, in the MFV MSSM, the mSUGRA scenario characterised by the 
input parameter values specified in Eq. (|2ip approximately maximises the cross section 
of process ([I]), which still comes out much below the threshold of observability at HERA 
and a future electron-proton supercollider. 

4.3 MSSM with non-minimal flavour violation 

Prior to presenting our NMFV-MSSM results, we explain our choice of input parame- 
ters. Scanning the eight-dimensional parameter space defined at the end of Section [3] we 
find that the following assignments, which we henceforth take as default, approximately 
maximise the cross section of process (JI]): 

tan/? = 8, M A o = 170 GeV, M = 475 GeV, M 2 = 655 GeV, 

M 3 = 195 GeV, A = 950 GeV, n = 345 GeV, A = 0.73. (22) 

These parameters are in accordance with the lower bounds on the squark masses of 
100 GeV and with the experimental lower bounds for the masses of the other SUSY 
particles [18) . They are also in accordance with the experimental bounds from b — *■ 57 
decay [18]. We calculate the MSSM spectrum using the program package FeynHiggs [22]. 

Figures [TTHH] all show the total cross section of process ([1]) in the NMFV MSSM. 
Specifically, Fig. [11] displays the dependence, also separately for the photoproduction 
and electroproduction contributions, while Figs. [T2"HT41 exhibit, for HERA experimental 
conditions, the dependencies on A, M , and M 3 , respectively. In Figs. d2] and [T3] also 
the contributions from loops involving gluinos, charginos, and neutralinos are shown sep- 
arately. 

From Fig. [11] we read off values of order 10~ 4 fb and 10 _1 fb for the cross sections at 
HERA and the future electron-proton supercollider mentioned above, respectively, which 
is discouraging in the case of HERA and challenging for the electron-proton supercollider, 
depending on its luminosity. 

From Fig. [12] we observe that the cross section strongly increases with A, by five 
orders of magnitude as A runs from to 0.73. For A > 0.73, the mass of the lightest 
squark is less than the lower limit of 100 GeV. The cross section is almost exhausted by 
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the gluino contribution. This may be understood by observing that the corresponding 
Feynman diagrams are enhanced by a factor of a s /a relative to those of the chargino and 
neutralino contributions. The neutralino contribution exhibits a A dependence similar to 
the full cross section, but is more than four orders of magnitude smaller. The chargino 
contribution oscillates about a mean value of 10~ 10 fb in the A range considered. For 
A <^0.3 (A ^0.3), it overshoots (undershoots) the neutralino contribution. 

From Fig. [13] we observe that the cross section strongly decreases with increasing value 
of M , by more than three orders of magnitude as M runs from 475 GeV to 2 TeV. The 
dominant role of the gluino contribution observed in Fig. [12] attenuates in the large- Md 
regime, where the chargino contribution gains influence. Nevertheless, there is a clear 
hierarchy among the gluino, chargino, and neutralino contributions for Mo ^ 700 GeV, 
the latter one being least important. For Mo ^700 GeV, the neutralino contribution 
exceeds the chargino one. The chargino and neutralino contributions exhibit minima at 
650 GeV and 850 GeV, respectively. 

From Fig.dHwe learn that the cross section decreases by about two orders of magnitude 
as M 3 runs from 195 GeV, the Tevatron search limit [TS], to 2 TeV. 

5 Conclusion 

In this paper, the photoproduction and photonic electroproduction of single top quarks in 
electron-proton scattering was analysed for the first time. The analysis was performed at 
one loop in the SM as well as the MSSM with minimal and non-minimal flavour mixing. 
In all three models, the cross section turned out to be too small to be measurable at 
HERA. The physics at HERA remains interesting because, in the light of our results, the 
single-top-quark-like events seen by the HI Collaboration might be a sign of physics not 
only beyond the SM, but also beyond the MSSM with conserved R parity. 
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Figure 1: Schematic representation of the hadronic process ([I]) explaining the four- 
momentum assignments. 
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Figure 2: Feynman diagrams of partonic subprocess (|2J) in the SM. Here, G is the charged 
Goldstone boson and di with i = 1, 2, 3 are the down-type quarks. 
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Figure 3: Additional Feynman diagrams of partonic subprocess (j2J) arising in the MFV 
MSSM. Here, df with i = 1,2,3 and s = 1,2 are the down-type squarks, and x% with 
i = 1,2 are the charginos. 





Figure 4: Additional Feynman diagrams of partonic subprocess arising in the NMFV 
MSSM. Here, g is the gluino, u a with a = 1, . . . , 6 are the up-type squarks, and x° with 
i = 1, . . . , 4 are the neutralinos. 
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Figure 6: Q 2 distribution of the cross section and its transversal and longitudinal parts 
in the SM under HERA experimental conditions. 
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Figure 7: Total cross section and its photoproduction and electroproduction parts in the 
MFV MSSM as functions of ^S. 
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Figure 8: Total cross section and its Higgs and chargino parts in the MFV MSSM as 
functions of tan (3 under HERA experimental conditions. 
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Figure 9: Total cross section and its Higgs and chargino parts in the MFV MSSM 
functions of m\/2 under HERA experimental conditions. 
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Figure 10: Total cross section and its Higgs and chargino parts in the MFV MSSM 
functions of A under HERA experimental conditions. 
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Figure 11: Total cross section and its photoproduction and electroproduction parts in the 
NMFV MSSM as functions of t/S. 
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Figure 12: Total cross section and its gluino, chargino, and neutralino parts in the NMFV 
MSSM as functions of A under HERA experimental conditions. 
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Figure 13: Total cross section and its gluino, chargino, and neutralino parts in the NMFV 
MSSM as functions of M under HERA experimental conditions. 
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Figure 14: Total cross section in the NMFV MSSM as function of M 3 under HERA 
experimental conditions. 



19 



